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ABSTRACT 


A  smooth  nonparametric  estimate  of  the  quantile  function  from  right- 

.1  i  a 

censored  data  is  given  by  Q^p)  =  h  Jg  Qn(  t)K(  ( t-p)/h)dt,  0  <  p  <  1,  where 
K  is  a  kernel  function,  h  is  the  "bandwidth,"  and  Qn  is  the  product-limit 
quantile  function.  This  report  describes  a  computation  procedure  for  data- 
based  selection  of  a  "best"  bandwidth  value  to  use  in  computing  Qn(p)  and 
for  obtaining  estimates  of  the  bias  and  standard  error  as  well  as  a  nonpara¬ 
metric  confidence  interval  for  the  true  quantile  value.  The  procedure  is  based 
on  the  bootstrap  method  for  right-censored  data.  A  listing  of  a  FORTRAN 
program  which  performs  the  necessary  calculations  is  provided,  and  examples 
of  the  procedure  are  given. 

Key  Words:  Nonparametric  quantile  estimation;  Right-censoring;  Percentile 
interval;  Bootstrap;  Bandwidth  selection. 


1.  INTRODUCTION 


Right-censored  data  arise  very  naturally  in  life  testing  and  reliability 
studies.  For  such  data,  it  is  important  to  be  able  to  obtain  good  nonpara- 
metric  estimates  of  various  characteristics  of  the  unknown  lifetime  distribu¬ 
tion.  This  report  concerns  the  computational  procedure  for  a  kernel-type 
nonparametric  estimator  of  the  quantile  function  of  the  lifetime  distribution 
from  right-censored  data.  This  estimator  was  suggested  by  Padgett  (1986), 
extending  the  complete  sample  results  of  Yang  (1985).  The  large  sample 
properties  of  the  estimator,  such  as  asymptotic  normality  and  mean  square 
convergence,  were  studied  by  Lio,  Padgett  and  Yu  (1986)  and  by  Lio  and  Padgett 
(1985). 

In  this  report,  a  procedure  for  calculation  of  the  kernel-type  quantile 
estimate  from  right-censored  data  is  described,  and  a  listing  of  a  computer 
program  in  FORTRAN  code  is  provided.  In  the  procedure,  the  bootstrap  is  used 
to  determine  the  approximate  "optimal"  bandwidth  values  to  use  for  the  kernel 
quantile  estimates  ("optimal"  in  the  sense  of  smallest  estimated  mean  squared 
error).  The  bootstrap  method  also  provides  estimates  of  the  bias,  mean  squared 
error, and  standard  deviation  of  the  estimator.  In  addition,  a  bootstrap  confi¬ 
dence  interval  for  the  quantile,  that  is,  a  percentile  interval,  is  obtained 
(see  Efron,  1982,  or  Efron  and  Tibshirani,  1986,  p.  68).  See  Padgett  and 
Thombs  (1986)  for  details  on  the  performance  of  this  computational  procedure. 

In  section  2,  the  estimator  is  given.  The  computation  procedure  using 
the  bootstrap  method  is  described  in  section  3.  Section  4  contains  the  listing 
of  the  FORTRAN  source  code  and  several  examples  on  the  utilization  of  the 


computer  program. 
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2.  THE  NONPARAMETRIC  QUANTILE  ESTIMATOR 


Let  X°,...,X°  denote  the  true  survival  times  of  n  items  or  individuals 
that  are  censored  on  the  right  by  a  sequence  U^,...,Un,  which  in  general  may 
be  either  constants  or  random  variables.  The  X?'s  are  nonnegative,  indepen¬ 
dent,  identically  distributed  random  variables  with  common  unknown  distribution 
function  Fq  and  unknown  quantile  function  Q°(p)  a  FQ^(p)  = 
inf{t:  Fq( t)  >  p},  0  <  p  <  1. 

The  observed  right-censored  data  are  denoted  by  the  pairs  (X^,A^), 
i-l,...,n,  where 

1  if  X°  <  U. 

i  -  i 

Xi  -  min{X°,U.},  &i  »  ■ 

.0  if  X?  >  Ui. 


Thus,  it  is  known  which  observations  are  times  of  failure  or  death  and  which 
ones  are  censored  or  loss  times..  The  nature  of  the  censoring  depends  on  the 


U^'s.  (i)  If  U^,...,Un  are  fixed  constants,  the  observations  are  time- 

truncated.  If  all  U^'s  are  equal  to  the  same  constant,  then  the  case  of 
Type  I  censoring  results,  (ii)  If  all  *  x°rj»  the  rth  order  statistic 
of  X°,...,X°,  then  the  situation  is  that  of  Type  II  censoring,  (iii)  If 

_ _ ,Un  constitute  a  random  sample  from  a  distribution  H  (usually  unknown) 

and  are  independent  of  X°,...,X°,  then  (X^,A^),  i=l,2,...,n,  is  called  a 
randomly  right-censored  sample.  The  observed  value  of  (X^,A^)  is  denoted  by 


A  popular  estimator  of  the  survival  function  l-FQ(t)  from  the  censored 
sample  (X^,A^),  i«l,...,n,  is  the  product-limit  estimator  of  Kaplan  and  Meier 
(1958).  The  product-limit  estimator,  which  was  shown  to  be  "self-consistent" 


by  Efron  (1967),  is  defined  as  follows.  Let  (Z^,A^),  i=l,...,n,  denote  the 
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ordered  X^'s  along  with  their  corresponding  A^s.  Then  the  product-limit 

estimator  of  1-F  (t)  is 
o 

1,  0  <  t  <  Zy 

"  'sSriA  zk-i  < 1  *  V k-2 . " 

1=1 

0,  zn  <  t. 

A  A 

The  product-limit  estimator  of  FQ(t)  is  denoted  by  Fn(t)  =  1-Pn(t),  and 

A  A 

the  size  of  the  jump  of  P  (or  F  )  at  Z.  is  denoted  by  s..  Note  that 

n  n  j  J  j 

Sj  -  0  if  and  only  if  Z_.  is  censored  for  j  <  n,  i.e.  if  and  only  if 

i 

A.  =  0.  Define  S.  =  I  s.  =  F  (Z.  .),  i=l,...,n-l,  and  S  *  1. 

J  1  j_l  J  “  1  +  1  n 

A  natural  estimator  of  Q°(p)  is  the  product-limit  (PL)  quantile 
function  Q^P)  *  inf{t:  Fn(t)  >  p}  (see  Sander,  1975,  or  Cheng,  1984,  for 
example).  Since  (^  is  a  step  function  with  jumps  at  the  uncensored 
observations,  it  is  desirable  to  obtain  a  smooth  estimator  of  Q°.  One  such 
estimator  is  the  kernel  estimator  defined  as  follows  (see  Padgett,  1986): 

Let  {h  a  hnJ  be  a  "bandwidth"  sequence  of  positive  numbers  so  that  hfi  -*  0 
as  n  ->  ®,  and  let  K  be  a  bounded  probability  density  function  which  is  zero 
outside  a  finite  interval  and  which  is  symmetric  about  zero.  Then  the  kernel 
quantile  function  estimator  is  given  by 


Qn(p)  -  h  1  jJ  Qn(t)K((t-p)/h)dt 
_1  n  S . 

*  h  1  l  Z.  J.1  K( (t-p)/h)dt,  0  <  p  <  1,  (2.1) 

i-1  1  5i-l 

where  Sq  s  0.  Here,  the  kernel  function  K  is  taken  to  be  the  triangular 

kernel  on  [-1,1],  i.e. 

'  l-|x|,  -1  <  x  <  1, 

K(x)  = 

v  0,  otherwise. 


This  function  satisfies  all  conditions  assumed  by  Padgett  (1986),  Lio, 
Padgett  and  Yu  (1986),  and  Lio  and  Padgett  (1985). 

Procedures  for  choosing  an  "optimal"  value  of  the  bandwidth  to  use  in 
calculating  (^(p)  from  a  given  right-censored  data  set  and  for  estimating 
the  standard  error  and  obtaining  an  approximate  confidence  interval  are  de¬ 
scribed  in  the  next  section.  These  computations  are  based  on  the  bootstrap 
for  right-censored  data,  and  a  FORTRAN  program  for  performing  them  is  listed 
in  section  4. 


3 .  COMPUTATION  PROCEDURES 

The  effective  performance  of  the  kernel  estimator  critically  depends  on 
the  choice  of  the  value  of  the  bandwidth  h,  which  is  a  "smoothing  parameter." 
If  h  is  too  small,  not  enough  smoothing  is  done,  and  the  estimate  will  be 
"rough,"  showing  features  which  do  not  represent  the  true  quantile  function. 

If  too  much  smoothing  is  done,  i.e.  h  is  too  large,  important  features  of  the 
true  function  may  not  be  evident.  Therefore,  for  a  given  set  of  right-censored 
data,  a  method  of  selecting  a  reasonable  value  of  h  to  use  in  calculating 
Qn(p)  at  each  desired  p  is  needed.  Ideally,  the  value  of  h  (which  will 
depend  on  p  in  general)  that  minimizes  the  mean  squared  error  of  Qn(p)  (or 
some  other  criterion)  should  be  used  in  computing  Qn(p)-  However,  that  value 
of  h  is  not  known  and  a  method  of  estimating  it  from  the  censored  data  is 
required.  Here,  the  bootstrap  method  of  estimating  the  mean  squared  error  of 
Qn(p)  as  a  function  of  h  is  employed.  Then  the  value,  h(p),  which 
minimizes  the  bootstrap  estimated  mean  squared  error  is  selected  as  the 
"optimal"  bandwidth  for  calculating  Qn(p).  The  estimate  of  the  mean  squared 
error  as  a  function  of  h  is  obtained  from  300  bootstrap  samples  taken  from 
the  given  right-censored  data  set  (x.,5.),  i=l,..., 


n.  This  "data-based 


method  of  choosing  the  bandwidth  value  was  mentioned  by  Padgett  (1986)  and  was 


further  investigated  by  Padgett  and  Thombs  (1986). 


The  estimate  of  Q  (p)  given  by  Qn(p)  in  equation  (2.1)  is  calculated 


using  the  estimated  bandwidth  value,  h(p),  obtained  from  the  procedure  above 
To  estimate  the  bias,  standard  error,  and  a  nonparametric  confidence  interval 


(i.e.,  a  percentile  interval),  1000  bootstrap  samples  are  generated  from  the 


given  right-censored  data.  For  the  jth  bootstrap  sample,  an  estimate  Q^..(p) 


is  obtained  from  formula  (2.1),  j=l,2, . . . ,1000.  The  bias  of  Qn(p)  is  then 


estimated  by 


BIAS  = 


1000 


1000 


Z  Qni(P)  -  Qn(p)  ' 

L  j=l  3  J 


and  the  standard  error  of  Q  (p)  is  estimated  from 

n 


SE 


/ 

1 

1000  2  (  i-1  Qnj(P)) 

\ 

999 

V 

(Qnj(p))  1000 

> 

Also,  an  estimate  of  the  mean  squared  error  of  Q^p)  is  obtained  from 


2  2 

MSE  -  SE  +  BIAS  .  To  calculate  a  nonparametric  confidence  interval  for 


Q  (p),  the  central  95-percentile  interval  is  obtained  by  selecting  the  25th 


value  of  the  ordered  bootstrap  estimates  Qnj(p)  as  the  lower  bound  and  the 


975th  value  in  this  ordering  as  the  upper  bound  (see  Efron,  1982,  or  Efron  and 
Tibshirani,  1986,  for  a  complete  discussion  of  this  type  of  interval). 


4.  THE  COMPUTER  PROGRAM  AND  EXAMPLES 
A  source  code  listing  of  the  FORTRAN  program  which  performs  the  computa¬ 


tions  described  in  section  3  is  given  in  Table  1.  That  is,  this  program, 
chooses,  based  on  the  minimum  bootstrap  mean  squared  error  criterion,  the 
bandwidth  value  to  use  for  calculation  of  Q  ( p ) ,  calculates  Q  (p),  calcu- 


lates  the  bootstrap  estimates  of  the  bias,  mean  squared  error,  and  standard 
error,  and  gives  the  approximate  95-percentile  interval  for  the  true  quantile. 
The  IMSL  (1982)  subroutine  GGUBS  for  generating  uniform  random  numbers  between 
zero  and  one  is  called  by  the  FORTRAN  program  and  must  be  available  for  use. 

Table  2  gives  an  example  of  input  data  for  the  program  in  Table  1.  The 
first  line  is  the  number  of  observations,  n  =  15.  Lines  2-16  give  the 
observed  (x^,S^)  pairs,  separated  by  spaces,  line  17  gives  the  number  k  of 
quantile  values  to  be  estimated  (1  <  k  <  10),  and  line  18  gives  the  k=4 
values  of  lOOp,  separated  by  one  or  more  spaces,  at  which  Qn(p)  is  to  be 
computed  (p  =  .05,  .10,  .25,  and  .50  in  this  example).  Table  3  shows  the 
output  for  the  data  in  Table  2.  This  output  requires  approximately  one  minute 
CPU  time  on  a  DEC  VAX  785  computer.  Tables  4  and  5  show  the  input  and  output 
for  the  mechanical  switch  data  of  Nair  (1984),  which  has  n=40  observations 
with  23  of  them  censored.  Quantiles  for  p=.05  and  p=.25  are  calculated. 
These  computations  require  about  3  minutes  CPU  time  on  a  VAX  785. 

Padgett  and  Thombs  (1986)  discuss  the  performance  of  this  procedure  of 
estimating  quantiles  from  right-censored  data.  This  nonparametric  procedure 
seems  to  perform  quite  well,  even  for  relatively  small  samples  of  size  n=20 


Table  1.  FORTRAN  Program  Listing  for  Quantile  Estimation 


C  ESTIMATION  OF  QUANTILE  USING  KERNEL  ESTIMATOR  FROM  RIGHT- 

C  CENSORED  DATA 

C  WITH  BOOTSTRAP  EST.  OF  BIAS,  VARIANCE,  AND  OPTIMAL  BANDWIDTH 

C  KERNEL  -  TRIANGULAR  ON  (-1,1) 

IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  Z( 300 ) , IDEL( 300 ) ,X( 300 ) ,S( 300 ) ,PF( 300 ) ,ZS( 300 ) 
1  IDS( 300 ) ,W( 300 ) ,P( 10 ) ,QNPP( 1000 ) 

REALM  W 
REALM  ISEED 
ISEED=22285.D0 
WRITE( 6 , 100 ) 

100  FORMAT ( '  QUANTILE  ESTIMATION  FROM  RIGHT-CENSORED  DATA  BY 
1  THE  KERNEL  METHOD  ' ) 

C  INPUT  SAMPLE  SIZE  -  INTEGER  <  300 
READ  * ,N 

C  INPUT  DATA  PAIRS  Z, DELTA  (SEPARATED  BY  SPACES) 

WRITE  (6,148)  N 

148  FORMAT ( '0SAMPLE  SIZE  =',  15/' 0CENSORED  SAMPLE:'/'  ',9X, 'X 
1  DELTA' ) 

DO  44  1=1, N 

READ  *,Z(I),IDEL(I) 

44  WRITE  (6,149)  Z(I),IDEL(I) 

149  FORMAT ( '  ',F15.6,I4) 

XN=N 

C  INPUT  NUMBER  OF  PERCENTILES  DESIRED  -  K 
READ  *,K 

C  INPUT  PERCENTILES  DESIRED  -  (UP  TO  10  VALUES  OF  P,  0<P<100) 
READ  *,(P(I) ,1-1, K) 

HNS= . 01 
HNINC=. 02 
HNMAX=  .75 

C  ORDER  Z'S  AND  DELTAS 
NMl-N-1 
DO  15  1=1, NM1 
IP1-I+1 
DO  15  J=IP1,N 
I F ( Z ( I) .LE.Z(J) )  GO  TO  15 
TEMP-Z(I) 

Z(I)-Z(J) 

Z  (  J  )  =TEMP 
ITEM=IDEL( I ) 

IDEL( I )=IDEL( J) 

IDEL- J )=ITEM 
15  CONTINUE 

DO  900  I 11=1 , K 
IFLAG=0 
NBSMAX-300 
PT?=?{ III )/100.D0 
QMSESAV =1 . D 1 5 
HN=HNS 


Table  1.  Continued 


1  DO  17  1=1, N 
ZS(I)-Z(I) 

17  IDS( I )=IDEL( I ) 

QBOOT=O.DO 

QBSQ=0.D0 

NBS=0 

14=1 

C  CALCULATE  PL  ESTIMATE  AND  S(I)=JUMP  SIZES 
3  DO  20  1=1, NMl 
PF( I )=1 .DO 
DO  20  J=1,I 

IF( IDS( J ) . EQ. 0 )  GO  TO  20 
XJ=J 

PF ( I ) =PF ( I ) * ( XN-XJ )/( XN-XJ+1 . DO ) 

20  CONTINUE 
PF(N)=0.D0 
S(1)=1.D0-PF(1) 

DO  21  1=2, N 

IMl=I-l 

21  S ( I ) =PF ( IMl ) -PF ( I ) 

IF(NBS.GT.O)  GO  TO  45 

C  CALCULATE  PL  QUANTILE  ESTIMATE 
SJ=0.D0 
DO  35  J=1,N 
SJP=SJ+S( J ) 

IF(S( J) .EQ.0.D0)  GO  TO  35 
I F ( ( PTP . LE. SJP ) .AND. ( PTP .GT.SJ ) )  GO  TO  36 
SJ=SJP 
GO  TO  35 

36  QNHATP=ZS ( J ) 

GO  TO  45 

35  CONTINUE 

C  CALCULATE  KERNEL  QUANTILE  ESTIMATE 

45  SJ=0.D0 

QP=0.D0 
DO  59  J=1,N 
A=  ( SJ-PTP ) /HN 
A=DMAXl ( A, -1 . DO ) 

SJ=SJ+S(J) 

XARG=( SJ-PTP )/HN 
B=DMINl ( XARG , 1 . DO ) 

IF( ( A.GE . 1 .DO ) .OR. ( B . LE.-l .DO ) )  GO  TO  59 
IF(B.LE.O.DO)  GO  TO  51 
IF(A.GE.O.DO)  GO  TO  52 
QP=QP+ZS ( J ) * ( B-A- ( B*B+A*A ) /2 . DO ) 

GO  TO  59 

51  QP=QP+ZS( J)*(B-A+(B*B-A*A)/2.D0) 

GO  TO  59 

52  QP=QP+ZS ( J ) * ( B-A- ( B*B-A*A)/2 . DO ) 

59  CONTINUE 
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Table  1 .  Continued 


IF(NBS.GT.O)  GO  TO  60 

CNP=<2P 

GO  TO  61 

60  QBOOT=QBOOT+QP 
QBSQ=QBSQ+QP*QP 

I F ( IFLAG.EQ.O)  GO  TO  61 

QNPP( 14 )=QP 

14=14+1 

61  NBS=NBS+1 

IF(NBS.GT.NBSMAX)  GO  TO  800 
C  GENERATE  N  UNIFORM  (0,1)  RANDOM  NUMBERS 
C  USING  IMSL  SUBROUTINE  GGUBS 
CALL  GGUBS(ISEED,N,W) 

C  GENERATE  BOOTSTRAP  SAMPLE 
DO  70  J=1,N 
WJ=DBLE(W(J) ) 

II=XN*WJ+1 . DO 
ZS( J)=Z( II ) 

70  IDS( J )=IDEL( II ) 

C  ORDER  ZS'S 

NM1=N-1 
DO  80  1=1, NM1 
IP1=I +1 
DO  80  J=IP1,N 

IF(ZS(I).EQ.ZS(J))  ZS( J)=ZS( J)+l.D-4 
IF(ZS(I).LT.ZS(J))  GO  TO  80 
TEMP=ZS(I) 

ZS(I)=ZS(J) 

ZS(J)=TEMP 
ITEM=IDS( I ) 

IDS( I )=IDS( J) 

IDS( J)=ITEM 
80  CONTINUE 
GO  TO  3 

800  XNS=NBS— 1 

QBIAS=QBOOT/XNS-QNHATP 

QMSE= ( QBSQ-QBOOT**2/XNS ) /( XNS-1 . DO ) +QBIAS**2 
IF(IFLAG.EQ.l)  GO  TO  860 
IF(QMSE.GE.QMSESAV)  GO  TO  850 
QMSESAV=QMSE 
HSAVE=HN 

I F ( I FLAG . EQ . 1 )  GO  TO  860 
850  HN-HN+HNINC 

IF(HN.LT.HNMAX)  GO  TO  1 
IFLAG=1 
NBSMAX-1000 
HN=HSAVE 
GO  TO  1 

860  QVAR=QMSE-QBIAS**2 


Continued 


WRITE( 6 , 103 ) 

103  FORMAT ( ' 0  P  QN(P)  BIAS  EST.  MSE  EST.  VAR  EST 

1  BANDWIDTH  USED' ) 

WRITE( 6 , 151 )  PTP,QNP,QBIAS ,QMSE,QVAR,HN 

151  FORMATS  ' ,F4.2,4D13.5,F9.2) 

SE=DSQRT ( QVAR ) 

WRITE( 6 , 152 )  SE 

152  FORMAT ( 'OESTIMATED  STANDARD  ERROR  OF  QN(P)  =',Dl2.5) 

DO  "/  51  1=1,999 
IPl-I+1 

DO  751  J=IPl,1000 

I F ( QNPP ( I ) . LE . QNPP ( J ) )  GO  TO  751 

TEMP=QNPP(I) 

QNPP ( I ) =QNPP ( J ) 

QNPP ( J ) =TEMP 
751  CONTINUE 

900  WRITE(6,153)  QNPP{25) ,QNFP(975) 

153  FORMAT ( '  APPROXIMATE  9 5-PERCENTI LE  INTERVAL: ' ,D13 . 5, 

1  '  TO  ' ,D13 . 5 ) 

STOP 

END 


Table  2. 


Example  Input  Data  for  Quantile  Estimation  Program 


15 

1.2837  0 

.6636  0 

.1827  0 

1.9805  0 

.1393  0 

.2796  1 

.6807  1 

.4247  1 

1.1301  0 

.3699  1 

1.9590  0 

.1404  0 

.1696  0 

.1912  0 
.4354  0 

4 

5  10  25  50 


n  «  sample  size 


?  pairs 


k  =  number  of  quantiles  to  estimate 
k  values  10 Op 
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Table  3.  Output  for  Data  in  Table  2 


QUANTILE  ESTIMATION  FROM  RIGHT-CENSORED  DATA  BY  THE  KERNEL  METHOD 
SAMPLE  SIZE  =  15 

CENSORED  SAMPLE: 

X  DELTA 

1.283700  0 

0.663600  0 

0.182700  0 

1.980500  0 

0.139300  0 

0.279600  1 

0.680700  1 

0.424700  1 

1.130100  0 

0.369900  1 

1.959000  0 

0.140400  0 

0.169600  0 

0.191200  0 

0.435400  0 

P  QN( P)  BIAS  EST.  MSE  EST.  VAR  EST.  BANDWIDTH  USED 

0.05  0.25144D+00  0.27371D-01  0.27924D-01  0.27174D-01  0.11 

ESTIMATED  STANDARD  ERROR  OF  QN(P)  =*  0.16485D+00 

APPROXIMATE  95-PERCENTILE  INTERVAL:  0.23801D+00  TO  0.57944D+00 

P  QN(P)  BIAS  EST.  MSE  EST.  VAR  EST.  BANDWIDTH  USED 

0.10  0.28883D+00  0.10934D+00  0.58516D-01  0.46560D-01  0.29 

ESTIMATED  STANDARD  ERROR  OF  QN(P)  =  0.21578D+00 

APPROXIMATE  95-PERCENTILE  INTERVAL:  0.22141EH-00  TO  0.10683D4-01 

P  QN(P)  BIAS  EST.  MSE  EST.  VAR  EST.  BANDWIDTH  USED 

0.25  0 . 77867D+00  0.37146D+00  0.20857EM-00  0.70595D-01  0.73 

ESTIMATED  STANDARD  ERROR  OF  QN(P)  »  0.26570D+00 

APPROXIMATE  95-PERCENTILE  INTERVAL:  0.36779D+00  TO  0.13626D+01 

P  QN( P)  BIAS  EST.  MSE  EST.  VAR  EST.  BANDWIDTH  USED 

0.50  0.14833D+01  -0.61632D+00  0.56409D+00  0.18424D+00  0.39 

ESTIMATED  STANDARD  ERROR  OF  QN(P)  =  0.42923D+00 

APPROXIMATE  95-PERCENTILE  INTERVAL:  0.56632D+00  TO  0.19792D+01 
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Table  4.  Input  Data  for  Mechanical  Switch  Failure  Example 


Cv> 

K\'.' 


t.' 


i 

I 

k: 


40 

1.151  0 
1.170  0 
1.248  0 
1.331  0 
1.381  0 
1.499  1 
1.508  0 
1.534  0 
1.577  0 
1.584  0 
1.667  1 
1.695  1 
1.710  1 
1.955  0 
1.965  1 
2.012  0 
2.051  0 
2.076  0 
2.109  1 
2.116  0 
2.119  0 
2.135  1 
2.197  1 
2.199  0 
2.227  1 
2.250  0 
2.254  1 
2.261  0 
2.349  0 
2.369  1 

2.547  1 

2.548  1 
2.738  0 
2.794  1 

2.883  0 

2.884  0 
2.910  1 
3.015  1 
3.017  1 
3.793  0 

2 

5  25 


S 

I 

Ifi 


Table  5.  Output  for  Mechanical  Switch  Failure  Data 


QUANTILE  ESTIMATION  FROM  RIGHT-CENSORED  DATA  BY  THE  KERNEL  METHOD 
SAMPLE  SIZE  -  40 

CENSORED  SAMPLE: 

X  DELTA 

1.151000  0 

1.170000  0 

1.248000  0 

1.331000  0 

1.381000  0 

1.499000  1 

1.508000  0 

1.534000  0 

1.577000  0 

1.584000  0 

1.667000  1 

1.695000  1 

1.710000  1 

1.955000  0 

1.965000  1 

2.012000  0 
2.051000  0 

2.076000  0 

2.109000  1 

2.116000  0 
2.119000  0 

2.135000  1 

2.197000  1 

2.199000  0 

2.227000  1 

2.250000  0 

2.254000  1 

2.261000  0 
2.349000  0 

2.369000  1 

2.547000  1 

2.548000  1 

2.738000  0 

2.794000  1 

2.883000  0 

2.884000  0 

2.910000  1 

3.015000  1 

3.017000  1 

3.793000  0 

P  QN(P)  BIAS  EST.  MSE  EST.  VAR  EST.  BANEWIDTH  USED 

0.05  0.16482D+01  0.43077D-02  0.12651D-01  0.12632D-01  0.05 

ESTIMATED  STANDARD  ERROR  OF  QN(P)  -  0.11239D+00 

APPROXIMATE  95-PERCENTILE  INTERVAL:  0.14995D+01  TO  0.19955EH-01 

P  QN(P)  BIAS  EST.  MSE  EST.  VAR  EST.  BANDWIDTH  USED 

0.25  0.21835D+01  -0.11022D-01  0.18869D-01  0.18747D-01  0.03 

ESTIMATED  STANDARD  ERROR  OF  QN(P)  -  0.13692D+00 

APPROXIMATE  95-PERCENTILE  INTERVAL:  0.18969D+01  TO  0.25470D+01 
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